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Abstract 

Among three typical energy scales, a neutrino mass scale {m v ~ 0.1 eV), a 
GUT scale (M GUT ~ 10 16 GeV), and a TeV-scale (M NP ~ 1 TeV), there 
is a fascinating relation of M NP ~ y/m u ■ M GUT . The TeV-scale, M NP , 
is a new physics scale beyond the standard model which is regarded as 
"supersymmetry" (SUSY) in this letter. We investigate phenomenology 
of SUSY SU(5) GUT with neutrinophilic Higgs, which realizes the above 
relation dynamically as well as the suitable magnitude of Dirac mass, m v , 
through a tiny vacuum expectation value of neutrinophilic Higgs. As a 
remarkable feature of this model, accurate gauge coupling unification can 
be achieved as keeping with a proton stability. We also evaluate flavor 
changing processes in quark/lepton sectors. 



PACS: 12.60.-i, 12.10.-g, 12.15.Ff 



1 Introduction 



There are three typical energy scales, a neutrino mass scale (m u ~ 0.1 eV), a GUT scale 
(Mqut ~ 10 16 GeV), and a TeV-scale (M NP ~ 1 TeV) which is a new physics scale beyond 
the standard model (SM) and regarded as supersymmetry (SUSY) in this letter. Among these 
three scales, we notice a fascinating relation, 

M 2 NP ~ m v ■ M GUT . (1.1) 

Is this relation an accident, or providing a clue to the underlying new physics ? We take a 
positive stance toward the latter possibility. 

As for a neutrino mass m„, its smallness is still a mystery, and it is one of the most im- 
portant clues to find new physics. Among a lot of possibilities, a neutrinophilic Higgs doublet 
model suggests an interesting explanation of the smallness by a tiny vacuum expectation value 
(VEV) pE]-[T5]. This VEV from a neutrinophilic Higgs doublet is of 0(0.1) eV which is the 
same as the neutrino mass, so that it suggests Dirac neutrino [31 HI El E] 3 Thus, the neutrino 
mass is much smaller than other fermions, since its origin is the tiny VEV from the different 
(neutrinophilic) Higgs doublet. Introduction of ^-symmetry distinguishes the neutrinophilic 
Higgs from the SM-like Higgs, where m v is surely generated only through the VEV of the 
neutrinophilic Higgs. The SUSY extension of the neutrinophilic Higgs doublet model is con- 
sidered in Refs. [7J [TTJ [121 US] • Since the neutrino Yukawa couplings are not necessarily tiny 
anymore, some related researches have been done, such as, collider phenomenology [81 fTU], low 
energy thermal leptogenesis [TTl [T2] . cosmological constraints |1 3^1 . and so on. 

On the other hand, SUSY is the most promising candidate of new physics beyond the SM 
because of a excellent success of gauge coupling unification. Thus, the SUSY SM well fits the 
GUT scenario as well as an existence of a dark matter candidate. 

There are some attempts that try to realize the relation in Eq. ( 11.11) . One example is to 
derive m u from a higher dimensional operator in the SUSY framework [16j. Another example 
is to take a setup of matter localization [T7] in a warped extra dimension [18]. These scenarios 
are interesting, but a model in this letter is much simpler and contains no additional scales 
other than M^p, m„, and Mqut- (For other related papers, see, for example, [T9l 120] .) 

In this paper, we investigate phenomenology of a SUSY SU(5) GUT with neutrinophilic 
Higgs (SU(B)h„) model proposed in Ref. [21]. E| Usually, SUSY neutrinophilic Higgs doublet 
models have tiny mass scale of soft ^-symmetry breaking (p,p' = O(10) eV in Refs. [11 [ [T2| 
IT5]). This additional tiny mass scale plays a crucial role of generating the tiny neutrino mass, 



* In Refs . [TJ [2l EJ [9j [TD1 [TTJ [T2] , Majorana neutrino scenario is considered through TeV-scale seesaw with a 
neutrinophilic Higgs VEV of 0(1) MeV. 

t A setup in Refs. [13] is different from usual neutrinophilic Higgs doublet models, since it includes a light 
Higgs particle. 

* A similar model was suggested in Ref. [22] , where lepton flavor violation was also roughly estimated. 
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however, its origin is completely unknown (assumption). In other words, the smallness of 
m v is just replaced by that of ^-symmetry breaking mass parameters, and this is not an 
essential explanation of tiny m„. This is a common serious problem exists in neutrinophilic 
Higgs doublet models in general. Notice, this problem can be solved by Ref. [2T] , where 
two scales of M GUT and M NP naturally induce the suitable magnitude of m u through the 
relation of Eq. (11.11) . and does not require any additional scales. The model contains a pair of 
new neutrinophilic Higgs doublets with GUT-scale masses, and the ^-symmetry is broken by 
TeV-scale dimensionful couplings of these new doublets to the ordinary SUSY Higgs doublets. 
Once the ordinary Higgs doublets obtain VEVs (v U)d ) by the usual electroweak symmetry 
breaking, they trigger VEVs for the neutrinophilic Higgs doublets of v Utd M NP / 'M GUT (~ m u ). 
Then, 0(1) Yukawa couplings of the neutrinophilic doublets to L N (L: lepton doublet, 
N: right-handed neutrino) give neutrino masses of the proper size. We can also obtain a 
GUT embedding of the SUSY neutrinophilic Higgs doublet model, which realizes the relation, 
m v ~ Vu^Mnp/Mqut, dynamically. As a remarkable feature of this model, accurate gauge 
couplings can be unified as keeping a proton stability. Flavor changing processes are also 
sensible aspect of this model. In general, flavor violation in charged lepton sector is related to 
that in quark sector because lepton doublet and right-handed down-type quark are contained 
in a same multiplet in SU(5) GUT. Particularly, neutrino oscillation directly contributes flavor 
violations in both sectors. It is one of our purposes to evaluate such flavor violating processes. 

This paper is organized as follows. In section 2, we review a SUSY SU(5)h„ model. In 
section 3 and 4, we discuss a gauge coupling unification, and investigate flavor violations in 
SUSY SU(5)h„ model. These sections are main parts of this paper. In section 5, we present 
a summary. 

2 SUSY 577(5) GUT with neutrinophilic Higgs 

Before showing a SUSY SU(5)h„ model [21], we show a SUSY neutrinophilic Higgs doublet 
model at first. This has a specific parameter region which is different from Refs. [TTl [T2l [T5] . 
We introduce Z 2 -parity, where only vector-like neutrinophilic Higgs doublets and right-handed 
neutrino have odd-charge. The superpotential of the Higgs sector is given by 

W h = fiH u H d + MH v H v i - pE u E v , - p'H v H d . (2.1) 

H v (H v i) is a neutrinophilic Higgs doublet, and H v has Yukawa interaction of LH U N, which 
induces a tiny Dirac neutrino mass through the tiny VEV of (H u ). This is the origin of 
smallness of neutrino mass, and this paper devotes a Dirac neutrino scenario, i.e., m„ ~ 
(H u ) = 0(0.1) eV. On the other hand, H v i does not couple with any matters. H u and H d are 
Higgs doublets in the minimal SUSY SM (MSSM), and quarks and charged lepton obtain their 
masses through (H u ) and (H d ). Note that this structure is guaranteed by the Z 2 -symmetry. 
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Differently from conventional neutrinophilic Higgs doublet models, we here take M as the 
GUT scale and p,p,p' 0(1) TeV. The soft ^-parity breaking parameters, p and p', might be 
induced from SUSY breaking effects (see below), and we regard p and p' as mass parameters 
of new physics scale, M^p = 0(1) TeV. Remind that usual SUSY neutrinophilic doublet 
models take p, p' = O(10) eV (for C(l) TeV E-terms) pU E21 E5J- This additional tiny mass 
scale plays a crucial role of generating the tiny neutrino mass however, its origin is just an 
assumption. Thus, the smallness of m v is just replaced by that of p and p' . This is a common 
serious problem exists in neutrinophilic Higgs doublet models in general. The present model 
solves this problem, in which two scales of Mqut and M^p induce the suitable magnitude of 
m u dynamically, and does not require any additional scales, such as 0(10) eV. This is one of 
the excellent points in this model. 

Amazingly, stationary conditions make the VEVs of neutrinophilic Higgs fields be 

PVu p'v d , ^ 

v v = ir , «„ = — . (2.2) 

It is worth noting that they are induced dynamically through the stationary conditions, and 
their magnitudes are surely of O(O.l) eV. Since the masses of neutrinophilic Higgs H v and 
H v i are super-heavy as the GUT scale, there are no other vacua (such as, v u ^ ~ v v y) except 
for v Uj d ^> v u y [T5] . Also, their heaviness guarantees the stability of the VEV hierarchy, 
v u> d ^ v v y, against radiative corrections [TU [15]. It is because, in the effective potential, H u 
and H y i inside loop-diagrams are suppressed by their GUT scale masses. 

The model suggested in Ref. [21] has the GUT scale mass of neutrinophilic Higgs doublets 
in Eq. fl2.ll) . so that it is naturally embeded into a GUT framework, and it is the SUSY SU (5)h„ 
model. A superpotential of a Higgs sector at the GUT scale is given by 

W|ut = M()trS 2 + AtrS 3 + HTi jj + ^ u ^ u - M X EE - M 2 $ v $ v , (2.3) 

where £ is an adjoint Higgs whose VEV reduces the GUT gauge symmetry into the SM. $ v 
(iv) is a neutrinophilic Higgs of ( ant i-) fundamental representation, which contains H u (H' v ) 
in the doublet component (while the triplet component is denoted as T„ (T v )). $„ and $ v 
are odd under the Z 2 -parity. H (H) is a Higgs of (anti-)fundamental representation, which 
contains H u (Hd) in the doublet component (while the triplet component is denoted as T (T)). 
The VEV of £ and Mo,i,2 are all of O(10 16 ) GeV, thus we encounter so-called triplet-doublet 
(TD) splitting problem. Some mechanisms have been suggested for a solution of TD splitting, 
but here we show a case that the TD splitting is realized just by a fine-tuning between (E) 
and Mi. That is, (E) — Mi induces GUT scale masses of T,T, while weak scale masses of 
H u , Hd- This is a serious fine-tuning, so that we can not expect a simultaneous fine-tuned 
cancellation between (E) and M 2 . Thus, we consider a case that the TD splitting only works 
in H and H, while not works in $^ and $„. This situation makes Eq. (12.31) become 

W e J f = pH u H d + MH u H u i + M'TT + M"T u f u . (2.4) 
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This is the effective superpotential of the Higgs sector below the GUT scale, and M, M', M" 
are of £>(10 16 ) GeV, while p = 0(1) TeV. 

Now let us consider an origin of soft ^-parity breaking terms, pH u H u > and p'H v H d in 
Eq.f l2.ip . They play a crucial role of generating the marvelous relation in Eq. ( 11.1 ft as well as a 
tiny Dirac neutrino mass. Since the values of p, p' are of order 0(1) TeV, they might be induced 
from the SUSY breaking effects. We can consider some possibilities for this mechanism. One 
example is to take a non-canonical Kahler of [S^(H u H v i + H v H d ) + h.c.]o. Where F-term of 
£ could induce the p- and p'-terms effectively through the SUSY breaking scale as in Giudice- 
Masiero mechanism |23j. There might be other models which induce the p- and p'-terms in 
Eq.( l2.ip except for introducing a singlet £. 

3 Gauge coupling unification and proton-decay 

In this section, we discuss a characteristic feature of the gauge coupling unification and the 
proton-decay in SUSY SU(5)h„ model by focusing on a role of T v and T v . As for the minimal 
SUSY £{7(5) GUT model, in order to unify the gauge couplings, mass of T and T should 
be lighter than the GUT scale as 3.5 x 10 14 GeV < W < 3.6 x 10 15 GeV due to threshold 
corrections [23] . However, to avoid the rapid proton decay, M' must be heavier than the GUT 
scale [M' > Mqut)- Hence, it is difficult to achieve both accurate gauge coupling unification 
and enough proton stability in the minimal SUSY £{7(5) GUT. 

A situation becomes different in the SUSY SU {5)h v model. In this model, a superpotential 
of the Yukawa sector is given by 

wy = \u^h + v^/dyM-# + f«Mj*» 

at the GUT scale, where % and j are family indices, ipi, (pi, and rji are 10-plet, 5-plet, and 
singlet in £{7(5) gauge group, respectively, which are written in terms of MSSM fields as 

A = {Q,, e-^Ui, (Vkm^Ej}, & = {(Vn^Dj, (Vd^Lj}, Vz = {e"^iV,}. (3.2) 

Since Yukawa couplings are written as 

f Uij = f„ <"" »,,■ ft, = (V^hfdM)^ Ui = l-/ 01 (3-3) 
the superpotential in this basis is given by 

VJy =f Ui QiUiH u + (V^ M )ijf dj QiDjH d + f di EiLiH d 

+ fu j (y KM )jiE i UjT — -fuiQiQiT + (V£ u )ij f dj UiDjT — (y^ M )ijf dj QiLjT 

- fv 1 {V D ) ij N i L j H u + f Vi (V D ) l3 N % D 3 T v , (3.4) 
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where CP phases, <p Ui and <f) Ui , are omitted, for simplicity. The terms from the fourth to 
seventh in Eq. (13.41) cause proton-decay, which also exist in the minimal SUSY SU(5) GUT. 
Thus, we should take M' > Mqut to avoid the rapid proton decay. Meanwhile the last term in 
Eq. (13.41) has nothing to do with the proton decay. Since T„ and T v contribute beta functions 
of SU(3) C x U(1)y, accurate gauge coupling unification is achieved with T v and T v threshold 
corrections with 3.5 x 10 14 GeV < M" < 3.6 x 10 15 GeV. Therefore, the SUSY SU(h) Hv model 
can realize not only the accurate gauge coupling unification but also the proton stability. 
Remembering that M is the GUT scale, 0(1) % tuning between M and M" is needed, but 
it can happen. Or, no tuning is required when one of couplings is of (9(0.01), for example a 
coupoing of S^H u H u i. 



4 Flavor changing processes 

Flavor changing in the lepton sector is related to that in the quark sector, since L and D 
are contained in a same multiplet in SU(5)h„- Where, mixing angles in Vp> are expected to 
be large, and masses of left-handed slepton and right-handed down-type squark get sizable 
radiative corrections in off-diagonal elements of flavor space. Leading log approximation makes 
the off-diagonal elements 

* -^fliV^kiiVn^Sml + Ao^og^, (4.1) 
( 5m ih * -g^/i(VS)H(Vb)jy(3mg + A,)log^, (4.2) 

where Mp is the Planck scale, mo and Aq are universal scalar mass and trilinear coupling in 
mSUGRA scenario. Equation (14. ip originates from a loop diagram of N and H u , where an 
energy scale in renormalization group equations runs from Mp to M (H v , H u mass). On the 
other hand, Eq. (14.21) is induced from a loop diagram of N and T u , which runs from Mp to 
M" (T u , T v mass). Notice that loop effects of Eqs. fl4.ll) and (14.21) are different from those in 
SU(5) with right-handed neutrinos (5 , f/(5)ijAr)§ In SU(5)rn model, neutrinos are Majorana 
and counterparts of Eqs. (14.11) and (14.21) are given by 

( 5m lh - -^^(^) fcl (^M)^(n./) im (^) mj (3mg + Ac) log (4.3) 

(SmDv ~ -±^f Uk (V*) kl (V D ) k] (3m 2 + A ) log^, (4.4) 

where M^ l is a diagonal Majorana mass of Ni (I = 1, 2, 3). A mass matrix of is diagonalized 
by a unitary matrix V M which does not appear in Eq. (l4.4p because M N[ is usually assumed 

§ See, for example, [25] 
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to be smaller than M' (> Mqut)- By comparing the SU(5)rn model with the SU(5)h v model, 
we can find an advantage point in the latter model. This is a predictability, that is, flavor 
changing processes are strongly predicted since there are no degrees of freedom of Vm- This 
means that flavor violations in charged lepton sector are directly related to those in the quark 
sector through the large flavor mixings in the neutrino sector. And, even if the mass matrix of 
the right-handed neutrinos is diagonal, magnitudes of <5m| and 5m 2 ^ in the SU (5)h v model are 
different from those in SU(5)rn model. For example, Eq.f l4.4p can be a few percent smaller 
than that in SU{5)h v model due to their log factors. The magnitude of the log factor is 
log -rrf 1z log H77, since M' must be larger than the GUT scale for the proton stability and M" 
must be of (9(10 14 ) GeV for the accurate gauge coupling unification. 

Let us show results of numerical analyses in n — > ej, r — > /17, and 6—7-57 processes. 
Figures [1] and [2] show correlations between branching ratios of B(6 — > S7) and B(/i — > cy) 
with tan/3 = 10 and A = 0. In FigJTJ m 1/2 is varied from 500 GeV to 800 GeV by 100 GeV. 
As for Figf2} m^j is varid from 500 GeV to 1000 GeV by 100 GeV. mo is varied from 200 
GeV to 1200 GeV by 100 GeV for each line. Here, the Higgs mass, calclated by FeynHiggs 
[271 EH |29l [30], is varied around 118 GeV which is not excluded by ATLAS [31] and CMS [32]. 
In Figs. [T]and El sin 2 26 l 1 3 is taken by and 0.01, respectively. We consider that spectrum of 
neutrinos is hierarchical, and the v T - Yukawa coupling is of 0(1). The current upper bound 
on B(/i — > e^y) is 2.4 x 10 -12 by the MEG experiment [33]. Figure [2] shows that large #13 is 
restricted in /i — > ej. In this parameter region, B(6 — > 37) does not change drastically because 
m^-dependence is larger than m - dependence. 

Figures [3] and [4] show correlations between B(6 — > 57) and B(r — > ^7), which parameters 
are the same as Figsfj] and El respectively. B(r — >• ^7) does not reach the experimental 
upper bound in this parameter region. (The experimental upper bound for B(r — > /17) is 
4.4 x 10~ 8 by BABAR experiment [33].) Note that a ratio of B(r — > /i7)/B(/i — > ej) depends 
largely on #13, where other neutrino oscillation parameters are fixed. When #13 becomes large, 
B(r — > /i7)/B(/x — > e'-f) closes to 10. This behavior is consistent with Ref. [22]. 

We do not consider r — > e'-f process because the experimental upper bound for B(r — > 67) 
is 3.3 x 10~ 8 [33] which is the same order as B(r — > ^7). Ratio of B(r — > ej)/B(r — > /17) is 
roughly proportional to (V / d)| 1 /(Vd)32 < 1- Hence, B(r — >■ fij) is more stringent constraint 
than B(r — > ej). 

Finally we comment on the Daya Bay experiment, which has measured a non-zero 6* 13 
[35]. The best-fit value is given by sin 2 26'i3 = 0.092, and such large mixing angle gives more 
stringent constraint in fi — > ej. Figure [5] shows m - and m!/2-dependence of B(/i — > ej) with 
sin 2 2^13 = 0.092. We can see that mi/2 should be larger than 2 TeV in order not to excess 
the experimental bound in Figj5j As for neutrinoless double beta decay, it is forbidden in our 
setup because neutrinos are Dirac fermion with lepton number conservation. 
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Figure 1: Contour plot of B(6 -> sj) and B(/x 
e7) with sin 2 26»i 3 = 0. Here, we take 200GcV < 
m < 1200GeV, 500GeV < m 1/2 < 800GeV, A = 
0, and tan/3 = 10. Experimental upper bound for 
B(^~> e-y) is 2.4 x 10~ 12 . 



Figure 2: Contour plot of B(6 37) and B(n 
e 7 ) with sin 2 26>i 3 = 0.01. We take 500GeV < 
m i/2 < lOOOGcV. Other parameters are the same 
as in Fig. Q] 




Bfr^prt 10' 11 10" 10 10" 9 10" 8 



Figure 3: Contour plot of B(6 ->• S7) and B(r 
/X7) with sin 2 26»i 3 = 0. Here, we take 200GeV < 
too < 1200GeV, 500GeV < m 1/2 < 800GcV, A = 
0, and tan/3 = 10. Experimental upper bound for 
B(r /i7) is 4.4 x 10~ 8 . 



Figure 4: Contour plot of B(b sj) and B(r 
M7) with sin 2 26>i3 = 0.01. We take 500GeV < 
m i/2 < lOOOGcV. Other parameters are the same 
as in Fig. [3] 
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Figure 5: B(£i — > e^) is plotted as a function of m 
with sin 2 26»i3 = 0.092. We take m 1/2 = 1 TeV, 1.5 
TeV, 2 TeV, and 2.5 TeV. 

5 Summary 

Among three typical energy scales, a neutrino mass scale, a GUT scale, and a TeV (SUSY)- 
scale, there is a marvelous relation of Eq. fll.lj) . In this paper, we have investigated phe- 
nomenology of a SUSY SU(5)h v model proposed in Ref. [21]. This model realizes the relation 
of Eq. fll.ip dynamically as well as the suitable m„ through a tiny VEV of neutrinophilic Higgs. 
At first, we have discussed the gauge coupling unification and the proton stability. Fascinat- 
ingly, the SU(h)n v can realize not only accurate gauge coupling unification but also enough 
proton stability simultaneously, which situation is hardly realized in usual four- dimensional 
SU(5) GUTs. Next, we have investigated correlations between b — >■ and /i — > ej, r — > fij. 
Notice that B{b — >• sj), B(fi — > ej) and B{r — > fi'j) are correlated directly through neutrino 
mixing in the SU(5)h„ model, which is an advantage of this model over the SU(5)rn model. 
As shown in Eq. (14.31) . additional unknown degrees of freedom, parameters in Vm, are needed 
in the latter model. Therefore, flavor changing processes are strongly predicted in the SU(5)h i/ 
model. As for the dependence of #13, £>(/i — > ej) depends largely on it, so that £>(/i — > ej) 
is strongly limited in large #13. On the other hand, we have shown that B{b — >■ S7) does not 
depend largely on # 13 . 
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